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ABSTRACT 



We analyzed the combined near-infrared spectrum of a star-planet system with thermal emission atmospheric models, based on the 
composition and physical parameters of the system. The main objective of this work is to obtain the inclination of the orbit, the mass 
of the exoplanet, and the planet-to-star flux ratio. We present the results of our routines on the planetary system HD 217107, which 
was observed with the high-resolution spectrograph Phoenix at 2.14 [im. We revisited and tuned a correlation method to directly 
search for the high-resolution signature of a known non-transiting extrasolar planet. We could not detect the planet with our current 
data, but we present sensitivity estimates of our method and the respective constraints on the planetary parameters. With a confidence 
level of 3-cr we constrain the HD 217107b planet-to-star flux ratio to be less than 5xl0~ 3 . We also carried out simulations on other 
planet candidates to assess the detectability limit of atmospheric water on realistically simulated data sets for this instrument, and 
we outline an optimized observational and selection strategy to increase future probabilities of success by considering the optimal 
observing conditions and the most suitable candidates. 

Key words. Planetary systems - Stars: individual: HD 217107 - Techniques: spectroscopic 
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1. Introduction 

The characterization of the over 500 detected exoplanets has 
now begun to take place. Most of the studies are carried out 
at optical and infrared wavelengths, because this is where the 
planetary reflected light and thermal e mission peak, respectively . 
The discovery of transiting planets (Charbonneau e t al.1 120001: 
iHenrv et"aT1l20 00) allowed astronomers to constrain new phys- 
ical parameters such as the radii and masses of the planets, 
which are not measurable by the radial velocity method alone. 
It is on these systems that in the last years the planetary atmo- 
sphere characterization has achieved the most exciting progress 
through the use of spectroscopy and broadband photometry with 
space telescopes. Examples are the identification of molecules 
such as water abs orption (e.g. Tinetti et al] 120071) or methane 
dSwain et al.l l2~008). or the o bservation of the th ermal emission 
. variation with orbital phase dKnutson et al.ll2007l) . 

Although great improvements in characterizing the compo- 
sition of transiting Hot-Jupiters have been achieved, they only 
represent about 20% of the known extrasolar planetsQ. The char- 
acterization of non transiting planets would require the direct 
detection of their light, but the very low flux ratios between 
the planets and their host stars makes a direct detection a very 
challenging goal. Secondary eclipse observations from Spitzer 
show that planet-to-star flux ratios can be as high as 2.5 xl0~ 3 
between 3.6 and 24 |im (e.g. iKnutson et alj[200a) . At 2.14 
the expected flux should be less than these values. Many authors 
have attempted a direct detection of the Doppler-shifted sig- 
nature in high-resolution spectroscopy from ground-based tele- 
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scopes. In the optical ICameron etal ] (l!999h fried to observe the 
starlight reflected from the giant exoplanet Tau-Bootisb, they 
found an upper limit to the albedo and radius using a least- 
squares dec onvolution method that is wel l described in the ap- 
pendices of ICollier Cameron et al. (2002), later the author re- 
peated the anal ysis on v An dromeda b (| Collier Cameron et al.l 
12002b . Recently iRodler etal .1 (120081 1201 Ol) searched in the visi- 
ble spectra of HD 75289Ab and Tau-Bootis b and found upper 
limits for their albedos using a model synthesis method. They 
constructed a model of the observation composed by a stellar 
template plus a shifted and scaled-down version of the stel- 
lar template to simulate the starlight reflected from the planet, 
these models were compared to the data by means of x 1 - In 
the near-infrared, several attempts have been made to detect 
Hot-Jupiters by trying to distinguish the planetary thermal emis- 
sion from the starlight ([Wiedemann et al. l l200Tt lL ucas & Roche 
2002; iBarnes et al.ll2007l 120081 1201 Oh . thev also found upper 
limits for the emitted flux of the planets. All these authors have 
used their own variation of a method based on the same princi- 
ple of separating the planetary and stellar spectra given their rel- 
ative Doppler shifts. Only recently. [Snellen et al.l(l2010l) claimed 
the detection of carbon monoxide from the transmission spec- 
trum of HD 209458 b during a ttansit observation by using high- 
resolution spectra; nonetheless, his technique required a transit- 
ing system. 

In this work we present an effort to constrain new physical 
parameters of the non-transiting Hot-Jupiter HD 217107 b. We 
attempt to trace its Doppler-shifted signature (estimated to be 
~10 -4 times dimmer than the star flux) with a correlation func- 
tion between high-resolution data and models of its atmospheric 
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spectrum. With positive detections this method would provide 
new information on its characteristics, such as its temperature, 
chemical composition, and the presence of chemical tracers as- 
sociated with life. At the same time, the method enables the cali- 
bration of high-resolution spectroscopic models for a larger sam- 
ple of planets that do not necessarily transit their parent star. 

In Section [2] we review the planetary system HD 217107; 
in Section [3] we describe the observations, data reduction, and 
calibration procedures; in Section |4]we detail the theoretical at- 
mospheric spectrum of the planet and the method used to extract 
and analyze the planetary signal and present the results of our 
data; in Section [5] we develop a strategy for the ideal data ac- 
quisition situation and simulate observations of other planetary 
systems; and in Section[6]we give the conclusions of our work. 



2. The planetary system HD 217107 

2.1. HD 21 7107b discovery 

HD 217107 is a main-sequence star that is similar to the Sun 
in mass, radius, and effective temperature; its spectral type, G8 
IV, indicates that it is sta rting to evolve into the red-giant phase 
dWittenmver et al l 120071). The presen ce of HD 217107b was 
first reported by iFischer et al.l (H999) through radial velocity 
measurements of the sta r, the detection was then confirmed by 
iNaef et all (l200ll) . Later. IFischer et al l (l200ll) id entified a trend 
in the residuals of the fit, and IVogt et al. (2005) postulated the 
existence of a third companion in an external orbit with a period 
of 8.6 ± 2.7 yr. The presence of this third o bject promoted the 
study of this system in subsequent surveys (iButler et al.ll2006t 
Wittenmv eret all 120071; IWright et al.l l2009). constraining more 
precisely the companions' orbital parameters. Table [1] summa- 
rizes the parameters used in this work. 



Table 1. Orbital parameters of HD 217107. 



Parameter 


Value 




References" 


Star: 








Spectral type 


G8IV 




W07 


7W(K) 


5 646± 


26 


W07 


K (mag) 


4.536 ± 


0.021 


C03 


d (pc) 


19.72 ± 


0.30 


P97 


M s (M Q ) 


1.02 ± 


0.05 


S04 


K s (ms" 1 ) 


140.6 ± 


0.7 


W07 


v g (kms -1 ) 


-14.0 ± 


0.6 


N04 


Planet: 








P (days) 


7. 12689 ± 


0.00005 


W07 


T p (3D) 


2 449 998.50 ± 


0.04 


W07 


e 


0.132± 


0.005 


W07 


m p sin i (M Jup ) 


1.33 ± 


0.05 


W07 


a (AU) 


0.074 ± 


0.001 


W07 


oj (deg) 


22.7 ± 


2.0 


W07 



07j; 

P97: |Perrvman & ESAldl997[) . S04: ISantos et all ( 120041) . 



2.2. Radial velocity 

The radial velocity of the planet, v p sin i, around the center of 
mass of the system is given by the reflex motion of the star: 




0.4 0.6 

Orbital Phase 

Fig. 1. Radial velocity curve of HD 217107 vs. orbita l phase . 
The crosses mark the observations of IWittenmver et all (120071) . 
which we used to compute this orbital solution. The boxes over 
the curve indicate the coverage of our observations, the filled 
boxes represent the runs utilized in the analysis, while the open 
boxes represent the discarded runs (details in Section H~2l . 



It depends on the mass of the star, m s ; the minimum mass 
of the planet, m p sin i; the projected radial velocity curve of the 
host star, v s (t) sin i (which in turn depends on the parameters 
T p , P, e, a), and K s ); and the inclination of the orbit, i, and also 
on the velocity of the center of mass of the system, v g , when mea- 
sured from Earth. Thus, the radial velocity curve of the planet 
is a distinctive curve in time, parameterized by the values sum- 
marized in Table [TJ where the only unknown parameter is the 
inclination of the orbit. Figure [TJ shows the radial velocity curve 
of the star owing to the interaction with HD 217107 b, phased 
over one orbit, with the origin in phase (<p — 0) at the time of pe- 
riastron. The radial velocity of the planet is proportional to this 
radial velocity curve (Equation [TJ. 



2.3. Flux estimate 

By simulating the spectra of the planet and its host star as black 
bodies, we can estimate the order of magnitude of the planet- 
to-star flux ratio as a function of wavelength. The black body 
emission, Fx(T), is determined by the surface temperature of the 
object. While for the star the temperature is well known from 
models (see Table [TJ, for the planet our best approximation is 
the equilibrium temperature 



T eq = 



1 - A 



1/4 



R s \ 1/2 , 



'eff- 



(2) 



For a reference value of the bond albedo of A = 0, we 
found an equilibrium temperature for HD 217107 b of T eq - 
1040 ± 19 K. Figure|2]shows the black body spectrum of the star 
and the planet assuming a radius between one and two Jupiter 
radii, which is the range of the radii for giant extrasolar planets 
measured to date. 

The planet-to-star flux ratio is given by 



Flux ratio = 



Fa(T planer) 
FA(T s tar) 



B A (T= 1040K) /fl ; , 
B A (T = 5646 K) \R S 



(3) 



v p (f) sin i = - v s (t) sin i 



x sin i. 



m„ sin i 



(1) 



At 2.14 |i.m, the flux ratio varies between 3x10 5 and 
1.5 x!0~ 4 from one to two Jupiter radii of the planet's radii, re- 
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Fig. 2. Black body emission of HD 217107 and HD 217107 b 
assuming a planet radius of 1.0 and 2.0 Jupiter radii. The verti- 
cal dashed line marks the waveband of our data (2.14 \im). The 
directly reflected light component has little contribution in the 
infrared and is thus omitted. 

spectively. For shorter wavelengths the flux ratio decreases, be- 
cause the star light dominates the emission spectrum. For longer 
wavelengths the net fluxes and thus the signal to noise ratio are 
lower. 

3. Observations and data reduction 

3.1. Observations 

We observed the planetary system HD 217107 in 11 nights 
between 200 7 Aug ust 14 and November 28 using Phoenix 
dHinkle et alJ l2003). a high-resolution near-infrared spectrome- 
ter at the Gemini South Observatory. 

The spectrograph has a 256 x 1024 InSb Aladdin II array 
with a resolving power of 10~ 5 |^m per pixel, the slit covers 14 
arc seconds in length. Its gain is 9.2 e~/ADU and it has a readout 
noise of 40 electrons. An argon hollow cathode wavelength cali- 
bration source is supplied with the instrument. Over 950 frames 
of the system were obtained in service mode, using the standard 
ABBA nodding sequences to easily remove sky emission, they 
cover a portion in the infrared spectral range from 2. 1 36 to 2 . 145 
|j.m (see Table [2]). We tuned the data acquisition after receiving 
the data from the first runs since the instrument was not fully 
characterized for use on the Gemini Telescope. For the first two 
nights, the exposure time was set to 25 seconds, whereas for the 
rest of the nights it was set to 80 seconds. We requested arc-lamp 
calibration exposures as well. 

3.2. Reduction 

We wrote our own interactive data language (IDL)0 routines 
for the data reduction and analysis, processing each night and 
slit position as an independent data set to minimize systematics 
caused by different atmospheric conditions or instrumental set- 
up. We used the flat-field images to identify hot pixels, marking 
a pixel as bad if it had a value beyond 3.5 sigma from the me- 
dian of the values of the nine subsequent pixels in its neighbor- 

2 http://www.ittvis.com/ProductServices/IDL.aspx 
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Table 2. Phoenix observations of HD 217107. 



Date 


Time on Target" 


Orbital Phase 


AV 


otatus 


UT 


min 




-1 

ms 




2007-08-14 


45 


0.30 


15.72 




2007-08-16 


45 


0.61 


12.29 




2007-08-22 


22 


0.43 


2.77 


rejected 


2007-08-26 


45 


0.99 


3.57 


rejected 


2007-10-02 


192 


0.16 


25.51 




2007-11-19 


96 


0.90 


2.61 


rejected 


2007-11-23 


96 


0.46 


2.79 


rejected 


2007-11-24 


96 


0.60 


13.90 




2007-11-25 


96 


0.74 


12.33 




2007-11-26 


96 


0.88 


3.84 


rejected 


2007-11-28 


96 


0.16 


10.23 





w Total exposure time of HD 217107. 
Radial velocity span of the star during the observing time. 
w See Sections RTTlandrOlfor details. 



hood. Bad pixels were masked in all further processing stages. 
Then, we divided the frames by a per-night master flat-field and 
subtracted their corresponding opposite A or B frame to remove 
bias and sky. Finally, we extracted the spectra from the frames 
with an IDL implementatio n^! of t he o ptimal spectrum extrac- 
tion algorithm described in Home (1986), this algorithm identi- 
fied cosmic ray hits, which were also masked from subsequent 
processing. 

3.3. Wavelength calibration 

First, we calibrated the wavelength dispersion using the ThAr 
lamps, identifying the line positions and stren gths in a high- 
resolution ThAr line atlas (lHinkleetal.1 12001). Because there 
was only one calibration lamp for each night, this solution rep- 
resented only a rough wavelength calibration, because there are 
(sub pixel) offsets in wavelength in the data. To reach the high 
precision needed for this work, we fine-tuned the calibration 
with a high-resolution spectrum of the Surfl to identify the tel- 
luric lines (identified as those present both in the solar spectrum 
and in an average spectrum of our data set). 

We constructed an average spectrum to increase the S/N ra- 
tio by aligning and adding the spectra of each night. To deter- 
mine the relative shifts, we selected within each set the first spec- 
trum as reference, while the rest were shifted (using spline inter- 
polations) to calculate the shift that minimized the root-mean- 
square of the correlation with the reference. The centers of fif- 
teen common absorption lines were identified in wavelength val- 
ues for the solar spectrum and in pixel position for our aver- 
age spectrum. The wavelength solution is obtained by fitting a 
second-order polynomial {A = co + c\ • p + c% • p 2 ) to the solar 
wavelength vs. the pixel position. Typical fitting coefficients are 
c = 2.145407, ci = -1.0305 x 10~ 5 , and c 2 = -1.8496 x lO" 10 . 
The dispersion of the residuals is RMS - 4.21 x 10~ 6 |J.m. 
No pattern is seen in the residuals. Pixels at wavelengths dom- 
inated by the identified telluric absorption lines were discarded 
from subsequent processing owing to their highly variable na- 
ture. About 65% of the pixels remained for the next analysis 
steps. 



3 http ://phy sic s . ucf . edu/ ~j h/ as t/sof t ware/optspecextr- 
0.3. 1/doc/index.html 

4 http://bass2000.obspm.fr/solar_spect.php 
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Fig. 3. Planetary spectrum blurring in the stellar template. Using 
our synthetic spectra of HD 217107 b we simulated the smear- 
ing of the planetary spectrum over one observing run. The dark- 
gray and light-gray lines denote the first and last spectra of 
HD 217107 b during an observing night. The relative shift ow- 
ing to the orbital motion of HD 217107 b is 40 kms 1 in this 
simulation. The bottom black line shows averaged the planetary 
spectra. 

4. Data analysis and results 

4.1. Correlation 

Because it is impossible to directly distinguish the planet's sig- 
nature from the stella r one i n a s ingle spectrum, fo l lowin g the 
idea of Demin g et alj d2000t) and Wiedem ann et alj d2001l) we 
searched for the planetary Doppler-shift signature through a cor- 
relation method between the (stellar-subtracted) residual data 
and a synthetic model of the planet's spectrum. To remove the 
stellar flux, we aligned the spectra for each set (Doppler-shifting 
them and using a spline interpolation) in a reference system in 
which the star remains at rest, and constructed a stellar tem- 
plate from the average of the set. Then, the stellar templates and 
the spectra are normalized dividing by their respective medians. 
Finally, the wavelengths of the stellar templates are shifted ac- 
cording to the orbital phase of the star in each individual spec- 
trum, and then the stellar template is subtracted from them. We 
avoided combining the different nights to obtain the stellar tem- 
plate, because it is highly probable that other systematics would 
be introduced. 

Because the planet is approximately a thousand times less 
massive than its host star, the planetary Doppler wobble is 
greater by the same order of magnitude (see Eq.[T], consequently 
the planetary signature will not be added coherently in the stellar 
template and thus appear blurred. The stellar template subtrac- 
tion leaves a residual spectrum that consists of the signature of 
the planet, which is slightly attenuated in the averaging process 
and immersed in Poisson noise. The blurring of the planet signa- 
ture (see Figure [3} is determined by the planetary velocity span, 
which in turn depends on the time span of an observation and 
the orbital phase at the time of the observation. Observations 
near inferior or superior conjunction provide the greatest radial 
velocity spans, while observations close to the greater elongation 
of the planet's orbit produce the smallest radial velocity spans, 
rendering the data useless. The rejected data sets in Table|2]were 
observed near greater elongation. 

For the high-resolution synthetic planetary spectra of 
HD 217107 b we used customized theoretical t hermal emis- 
sion models of its atmosphere (model described in lFortnev et al.l 



l2005Ll2006ll2008h at three different distances from the star to ac- 
count for the non-negligible eccentricity of the planet. The mod- 
els are cloud-free, with solar metallicity, gravity g = 20 ms~ 2 , 
and the molecular abundances are those appropriate for chemical 
equilibrium. At these effective temperatures, the main absorbing 
molecules are H2 O, CH4, CO, and CC>2- Th e c hemistry is de- 
scribed in detail in iLodd ers & Feglev (2002) and Visscher et al.l 
(120061) . We empirically characterized the instrumental resolu- 
tion through the analysis of the emission lines in the calibration 
lamps. We convolved the model spectra by the instrumental res- 
olution, which we determined to be A/AA w 40 000. Then, for an 
assumed value of sin i, the synthetic spectra are Doppler-shifted 
to mimic the radial velocity of the planet at the time that the data 
frame was obtained. The correlation degree, C(i), is calculated 
according to the formula 

m Ly z;.K-^)-M-^o) (4) 

In this equation we used the notation fy for the value of the 
function at the pixel j of the spectrum k, and for the mean 
value of the function in the spectrum k. Here, "r" refers to the 
residual spectrum while "t(/)" to the shifted planetary model 
spectrum, with Nk the number of pixels in spectrum k and N the 
total number of spectra. The denominator in the expression nor- 
malizes the correlation, and thus a value of 1.0 would indicate a 
perfect correlation. 

We thus produce a curve of the correlation degree vs. the in- 
clination of the orbit, evaluated in the range < i < n/2. A 
positive value of this function indicates that the data spectrum 
resembles that of the model, while a negative one suggests anti- 
correlation. As consequence of the random nature of the Poisson 
noise, the value of the correlation between the residual spectra 
and the models should be close to zero, except when the adopted 
i matches that of the planetary system. Therefore, an apprecia- 
ble peak in the correlation curve would represent a successful 
detection of the planetary signature and immediately indicates 
the value of i. By constraining the inclination with this method, 
the mass of the planet would be immediately determined via Eq. 
1. 

4.2. Data results 

For this analysis, we excluded the nights where the velocity span 
of the star was less than 10 ms" 1 (column 4 of Table |2]i since 
they do not represent any significant improvement in the results, 
because the shift of the planet (~ 8.3 kms~') is not significantly 
higher than the instrumental resolution (~ 7.5 kms -1 ). Figure|4] 
(Top panel) shows the correlation curve derived from our data as 
a function of sin i. The degree of correlation found was close to 
zero at all inclinations, and we do not distinguish any identifiable 
positive peak that could indicate an atmosphere with absorption 
features resembling those of the models. 

4.3. Planet-to-star flux ratio fitting 

While the inclination determines the maximum of the correlation 
curve, the planet-to-star flux ratio (F p /F s ) is the main physical 
parameter bounded to the magnitude of the correlation. In this 
section we determine the most probable values in the parameter 
space [F p /F s , sin i], which gave rise to our result, and estimate 
the statistical significance of the value of the correlation reached. 
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We searched for the best fitting values comparing our data results 
(Fig. 3] Top) with "synthetic" correlation curves. We generated 
the synthetic correlation curves by recreating our observations, 
adding a synthetic planetary spectrum, with known inclination 
and planet-to-star flux ratio, according to the following scheme: 
Step 1 : We rearranged the order of the data set with random per- 
mutations within each night, but kept the original order of the 
dates of the observations. As a consequence, any real planet sig- 
nature disappeared, but the noise level of the data was conserved. 
Step 2: Using the atmospheric models of the planet, we injected 
a synthetic spectrum in the scrambled data set, Doppler-shifted 
and with a relative flux according to specific values of sin i and 
F p /F s , respectively. For simplicity, we adopted a constant F p /F s 
along the orbit. 

Step 3: We processed these synthetic data through the same rou- 
tines as in our original data (section |4~TV We then iterated for a 
grid of values in the ranges: < i < n/2 and 10~ 5 < F p /F s < 
10~ 2 , obtaining a set of synthetic correlations for sin i and F p /F s . 

Once we obtained these models, we searched for the best-fit 
parameters through a x 2 minimization between the data corre- 
lation curve and the synthetic correlation curves, generating a 
goodness-of-fit map (Fig. |4] bottom panel). 

In addition, we used a bootstrap procedure to calcu- 
late false-alarm-probability limits for this map. Following 
ICollier Cameron et al.l (120021) . we determined the frequency 
with which the correlation degree exceeds a given value as a 
result of noise in the absence of a planet signal. The routine con- 
sists of performing a random permutation of the data sets and the 
subsequent data analysis (steps 1 and 3 of previous paragraphs) 
which we repeated a large number of times (~5 000), record- 
ing the correlation curve after each trial. This set of correlation 
curves represents the correlation found in the absence of a plan- 
etary signal, and, because it is created from the data themselves, 
defines an empirical probability distribution that includes both 
the photon statistics and instrumental systematics. 

Then, at each inclination, we stacked and sorted the values 
of the correlation in increasing order. We determined the 1, 2, 3, 
and 4-cr false-alarm confidence levels as the value of the corre- 
lation degree at the 65, 90, 99 and 99.9 percentiles of the trials. 
They represent the signal strengths at which spurious detections 
occur with 35, 10, 1, and 0.1 percent false-alarm probability re- 
spectively, at each value of the inclination. This allows us to as- 
sess the probability of obtaining a certain correlation degree in 
the absence of planetary emission. 

Fig. |4] (Bottom panel) shows the probability map for 
HD 217107b. The best fit occurs at sin / = 0.838 and F p /F s = 
3.6 X 10~ 3 , although the relative improvement in^- 2 against the 
surrounding parameters is shallow. This value disagrees with the 
maximum value of the correlation curve (Fig. |4]Top), and fur- 
thermore, the bootstrap results indicate that this value is below 
the 3-cr confidence limit of the signal not being a false positive. 
Also, this F p /F s is much higher than the predicted value from 
Sec. 1231 (between 3xl0~ 5 and 1.5xl0~ 4 ). The disagreement of 
the results of the top and bottom panel in Fig.|4]suggests that sys- 
tematics remain after the data reduction, while the strength of the 
correlation value, two orders of magnitude above the expected 
flux ratio, indicates that this result is not realistic. The 3-cr con- 
fidence limit only allows us to establish an upper limit in the 
flux ratio at 4-5 xlO -3 for inclinations greater than sin i = 0.6. 
In conclusion, we cannot state the detection of HD 217107 b. 
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Fig. 4. Top: Correlation result for our data set as a function 
of sin/. The correlation remains flat along every inclination 
without any distinctive peak, the maximum value is reached at 
sin / = 0.71. Bottom: Goodness of fit, ;^ 2 -map, of the correlation 
models to the data. The horizontal and vertical axes refer to the 
fitting parameters sin i and F p /F s respectively, at which the syn- 
thetic planetary spectrum was added in the correlation models 
creation, from which we calculated the minimum squares 0f 2 fr ). 
We plot x 2 relative to the best fit (Ay 2 fr = x 2 fr ~ X^) using the 
function exp(-a' • Ax 2 ). The gray scale denotes the goodness of 
fit, from black for the best fit (at xtm)' t0 white f° r the poorest 
fit. The plotting parameter, a, just enables a good contrast in the 
plots (the same value was used for all plots). Additionally, we 
determined with bootstrap procedures the solid lines (bottom to 
top) that mark the (1, 2, 3 and 4-cr) levels of false-alarm prob- 
ability. The white cross marks the best fit at sin i = 0.84 and 
Fp/Fs = 3.6 x 10 -3 , situated below the 3-cr confidence level. 

5. Future prospects 

5.1. Observational strategy 

Although our current data do not enable us to claim the detec- 
tion of HD 217107 b, we identified a strategy to maximize the 
chances of a successful detection. This involves selecting suit- 
able candidate systems and precisely choosing the phasing and 
span of the observations. To exemplify the advantages, we sim- 
ulated realistic observations of other planetary systems. 

We limited our sample to the currently known extrasolar 
planets without transits^ observables from Gemini South, for 
the first semester of the year. Even though we constrained our- 
selves to the instrument we characterized and to a fixed time 
span, the purpose of these simulations is to provide one success- 
ful detection with our method. It is plausible that considering the 
full extent of possibilities, stronger signals can be acquired. The 
improvement in a detection, limited by purely photon noise, can 
be quantified by the planet-to-star flux ratio and the stellar flux, 
according to the expression (fluxes in number of photons) 

(F p /F s ) ■ F s 

F /No i S e = pl * {Fp/Fs) ■ <F S . (5) 

Then, for example, the spectrometer CRIRES with four times the 
wavelength coverage of that of Phoenix, has twice the sensitiv- 

5 http://exoplanet.eu 
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Fig. 5. Correlation curves and ^ 2 -maps of synthetic data of HD 179949. A synthetic planetary signal was injected in the spectra 
with parameters: sin i = 0.77 and F p /F s = 0.003. Left: Using our observing strategy. Right: random distributed observing dates. 
The routine successfully recovers the signal at sin i = 0.78 and F p /F s = 2.8 x 10 3 in both cases, although, when using our strategy, 
the correlation degree is stronger, and the parameters are better determined compared with the right panel. 



ity of Phoenix. We decided to simulate the Phoenix instrument, 
since it is well characterized by our group, while other instru- 
ments should present their own systematics, which are hard to 
quantify. 

We simulated the planets as if the strength of the high- 
resolution absorption features were the same as that of our mod- 
els for HD 217107b, but with the corresponding F p /F s (esti- 
mated as in Sec. 12.31 ). A caveat for this assumption is that the 
strength of the l ines is not very clear i n pla nets that exhibit ther- 
mal inversions. iBurrows et al.l (120081) and iFortnev et al.l (|2008) 
suggest that the emission features should be weaker. 

The target selection criteria are based first on the radial ve- 
locity span of the planet, where we set a lower limit cutoff of 7.5 
km s' 1 (equivalent to the FWHM of the instrument spectral reso- 
lution) for a three-hour observing run if the orbit was at i = 30°. 
Second, we look for higher apparent brightnesses of the stars for 
better signal-to-noise ratios. Table|3]lists two of the better suited 
selected planetary systems (HD 217107 listed for comparison). 
A brighter K-band magnitude of the star improves the signal-to- 
noise ratio, while a smaller semi-major axis involves a higher 
radial velocity span, which enables a greater Doppler shift of the 
planet spectra during the runs and at the same time favors higher 
planet-to-star flux ratios. 

To simulate the observations, we recreated the same instru- 
mental settings of our data, but carefully selected the observ- 
ing schedule. For each one of the nights in the period and re- 
stricted to air masses under 1.5, we selected the three-hour range 
that gives the maximum velocity span. We recorded then, the 
radial velocity spans for each night, and chose those with the 

Table 3. Favorable targets for Gemini South. 



Target 


a 


M K 


Kp 




AU 




km s" 1 


HD 179949 


0.045 


4.94 


158.23 


Tau Boo 


0.046 


3.51 


150.62 


HD 217107 


0.073 


4.54 


112.28 



Notes. Planet's radial velocity amplitude for an orbit with sin i = 1. 



biggest spans. We used the solar spectrum to simulate the stellar 
spectrum, while for the planetary component we used the atmo- 
spheric models of HD 2 1 7 1 07 b added with a given planet-to-star 
flux ratio and inclination. Each component is Doppler-shifted 
according to the orbital parameters. Finally we added Poisson 
noise to the spectra, according to the signal-to-noise correspond- 
ing to the magnitude of the target. The synthetic data were pro- 
cessed in exactly the same way as our original data. 

5.2. Simulations 

In our first test, we present two simulations of an observing cam- 
paign on a target with the physical parameters of HD 179949 
to show the improvements of our observing strategy in contrast 
with a regular observation. The given parameters are sin i = 0.77 
and F p /F s = 3 x 10~ 3 . Figure left shows the simulation fol- 
lowing our observing strategy, while Fig.|5]right shows the sim- 
ulation selecting random observing dates. In both cases the cor- 
relation curves (top panels) mark the inclination of the synthetic 
orbit with a increment in the correlation degree near sin i = 0.77, 
while the ^ 2 -maps (bottom panels) effectively indicate the best 
fit at sin i = 0.78 and F p /F s = 2.8 x 10~ 3 . 

We identify the main differences between these two simula- 
tions: First, given the larger radial velocity spans when imple- 
menting our strategy, the planetary spectrum is more blurred in 
the stellar template and consequently less reduced when the tem- 
plate is subtracted, the planetary spectrum signal is thus stronger 
in the residual spectrum, which increases the correlation degree. 
As consequence of these greater correlation degrees, all confi- 
dence levels are generally lower, because it is less probable to 
reach this correlation degree by chance in the no-planet case, and 
lastly, the^ 2 -map peak is much better determined. The improve- 
ment is reflected more in the distinction of the best fit against 
other values of the parameter space than in the distinction against 
the no-planet case. 

In another simulation, we recreated the planetary system 
Tau Boo as close as possible to its real physical characteris- 
tics (Figure |6j. Tau Boob had the parameters sin / = 0.82 and 
F p /F s = 4 x 10~ 4 , our routines returned the best fit: sin i = 0.79 
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sin(z) 



Fig. 6. Similar to Fig. [4] correlation curve (top panel) and 
^- 2 -map (bottom panel) of a simulation of the planetary sys- 
tem Tau Boo, with an injected companion at sin i = 0.82 and 
F p /F s = 4 x 1(T 4 . 

and FplF s = 3.6 x 1CT 4 , slightly underestimating the values. 
Nevertheless, the ^ 2 -map shows an improvement in the region 
near the injected inclination and flux ratio. The bootstrap re- 
sults set the 3-cr confidence limit near F p /F s 1.5 xlCT 4 (for 
sin i > 0.5), indicating a detection with 99% confidence. 



recovered the planetary signature in both cases, but showing an 
improvement in the correlation degree, precision in the ^ 2 -map, 
and lower cr limits when using our observing schedule. Finally, 
we performed a realistic simulation of the planetary system Tau 
Boo, and successfully detected its signature. 

In conclusion high-resolution instruments like Phoenix are 
capable of detecting extrasolar planet Doppler-shifted signals 
with flux ratios as low as 10 4 with this method if we perform 
a careful treatment of the systematics (approaching the photon 
noise limit), if we count with appropriate theoretical models, 
and if we follow an optimized scheme in the data acquisition. 
Furthermore, using other instruments like CRIRES or NIRSPEC 
could increase the confidence of the results. Since our simula- 
tions exclude systematics effects specific to the instrument, an 
adequate treatment to remove them would be necessary. Our op- 
timal observing strategy tends to select observations at superior 
conjunctions of the planet's orbit, capturing the highest amount 
of light possible from the planet and at the same time cover- 
ing the highest radial velocity span for a determined time extent. 
Refinements of this technique will involve the optimization of 
the distribution of time designated to the length of an observing 
run vs. the number of nights of observation, while adding phase- 
dependent functions of the planet's brightness to account for the 
changing observed portion of the day/night side of the planet, 
and for different amounts of irradiation in eccentric orbits, will 
increase the accurateness of the fitted parameters. 

Acknowledgements. Patricio Cubillos and Patricio Rojo are supported by the 
FONDAP Center for Astrophysics 15010003, the center of excellence in 
Astrophysics and Associated Technologies (PFB06) and the FONDECYT 
project 11080271. 



6. Discussion and conclusions 

Because the instrument was not well characterized at the time 
and our service-mode observational strategy had to be adapted 
after the first few observing windows, the data for HD 217107 
were not as sensitive as expected. The correlation curve was fea- 
tureless for all inclinations and with values close to zero, with a 
maximum at sin/ = 0.71. By fitting the sine of the inclination 
and the planet-to-star flux ratio through least-squares, we found 
the best-fit parameters of sin i = 0.84 and F p /F s = 3.6 x 10~ 3 at 
a level below our 3-cr confident limit. As a consequence of the 
faint features in the results, the disagreement between the peak in 
the correlation (Fig. |4]Top panel) and the most probable value of 
sin i (Fig.|4]Bottom panel), and the higher than predicted F p /F s , 
we could not claim a detection of HD 217107 b with our current 
data. Given the results of the bootstrap procedure, we reject the 
flux ratio of HD 2 1 7 1 07 b to its host star to be over 5 x 1 3 (3-cr 
confidence). We attribute these results to the absence of an ideal 
strategy in the data acquisition at the time of the observations 
and a needed further treatment of the instrument systematics. 

We could not detect HD 217107 b, but defined the outlines of 
future campaigns by carefully defining a candidate selection cri- 
terion and an observational strategy. We conclude that the best- 
suited candidates for this technique are those in very close or- 
bits, which allow the planets to have high orbital velocities and 
higher planet-to-star flux ratios. We propose an observing strat- 
egy where for the period of observations we specifically select 
the nights with maximum radial velocity spans. To explore the 
capabilities of our routines, we simulated other planetary sys- 
tems as observed by the Phoenix spectrograph, with the same 
number of hours and an appropriate schedule of observations . 
The system HD 179949 was recreated, contrasting the use of 
our observing strategy with a regular observation schedule. We 
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